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The  structure  of  spontaneous  decay  of  atoms  in  the  vicinity  of  a  surface  is 
shown  to  be  determined  by  spatial  symmetries.  The  spontaneous -decay  operator 
for  a  degenerate  two- level  atom  is  derived,  and  with  symmetry  considerations  the 
nvimber  of  free  parameters  is  reduced  to  two.  Only  the  dimensionless  and 
normalized  inverse  lifetimes  bj|  and  bj^  for  a  parallel  and  perpendicular  dipole 
moment  with  respect  to  the  surface,  respectively,  enter  the  expression  for  the 
relaxation  operator  for  any  atom  near  any  surface.  These  two  parameters 
incorporate  the  atom- surface  distance  dependence  of  all  Einstein  coefficients 
for  spontaneous  decay,  and  all  optical  properties  of  the  substrate  material.  It 
is  shown  that  the  specific  features  of  spontaneous  decay  are  mainly  geometrical, 
and  a  consequence  of  symmetries  of  the  vacuum  radiation  field,  irrespective  of 
the  presence  of  the  atom.  With  an  example  it  is  shown  how  the  parameters  b||  and 
bj^  can  be  calculated  in  a  particular  case. 


An  excited  atom  will  decay  to  lower  states,  until  it  reaches  its  ground 
state.  This  process  of  spontaneous  decay  is  accompanied  by  the  emission  of 
fluorescent  photons  and  the  loss  of  internal  energy  of  the  atom  equals  the 
energy  gain  of  the  radiation  field.  Therefore,  the  processes  of  spontaneous 
decay  and  spontaneous  emission  of  radiation  are  related  through  energy 
conservation.  We  shall  consider  a  dipole -allowed  transition  between  a 


degenerate  excited  level  e  and  a  ground  level  g  (also  possibly  degenerate) ,  and 

indicate  the  atomic  wave  functions  of  the  multiplets  by  |j  m  >  and  |j  m  >, 

e  e  6  6 

respectively.  The  upper  state  has  2jg  +  1  magnetic  substates,  which  are  coupled 


by  the  atomic  dipole  moment  operator  ^  to  the  2j  +  1  ground  states.  If  we 

O 

denote  by  the  energy  separation  between  the  levels,  then  the  expression  for 


the  Einstein  coefficient  for  spontaneous  decay  from  e  to  g  reads 


A 


f 


o' 


(1.1) 


in  terms  of  the  reduced  matrix  element  of  the  dipole  operator.  An  atom  in  state 
|j^m^>  shall  decay  exponentially,  with  a  lifetime  1/A^,  to  the  various  ground 
states  at  a  rate  A^  times  the  population  of  after  completion  of  the 

decay  the  emitted  radiation  energy  equals  An  important  feature  of  this 

process  is  that  the  relaxation  constant  for  the  decay  of  is  independent 

of  the  magnetic  quantum  number  ra^.  That  this  must  be  so  is  a  consequence  of 
symmetry,  as  can  be  understood  as  follows.  If  we  rotate  an  atomic  wave  function 
|j^m^>  over  the  Euler  angles  a,/3,7,  then  the  transformed  wave  function  can  be 
expressed  as  [1] 


\tMi' 


I  \  I 

f^(a,fi.y)\i^m^>  -  )  |j^m;>  (a, ^.7)  . 

e  e 


(1.2) 


m' 
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where  Pg^(a,^,7)  indicates  the  rotation  operator,  and  D 


(jJ 


is  the  rotation 

matrix.  Spontaneous  decay  of  spontaneous  emission  of  photons,  is 

brought  about  by  a  coupling  of  the  atomic  dipole  moment  to  the  electromagnetic 
field.  But  since  the  electromagnetic  vacuum  (empty  space)  is  isotropic,  a 
rotate  state  P  (a,^,7)|j  m  >  must  decay  in  the  same  way  as  the  original  state 
|i  m  >,  for  all  rotation  angles  a, ^,7.  From  (1,2)  it  follows  that 

G  G 

^m^>  is  a  superposition  of  all  states  therefore  this  can 

only  hold  if  the  relaxation  constant  for  |jg™g>  independent  of  m^.  From  a 

different  point  of  view,  we  can  say  that  the  quantum  number  m^  refers  to  a 

particular  choice  of  the  quantization  z-axis.  Rotating  an  atomic  state  |jg™g> 
is  then  equivalent  to  a  change  of  quantization  axis.  In  isotropic  space  the 
choice  of  the  z-axis  has  no  significance,  and  consequently  the  decay  rate  of 
(j^m^>  must  be  independent  of  this  choice,  and  therefore  independent  of  m^. 

Let  us  now  consider  an  atom  which  is  positioned  near  a  surface.  We  choose 

the  z-axis  (arbitrarily)  as  the  normal  to  the  surface.  The  region  z  <  0  is 

filled  with  an  optically- reflecting  material,  like  a  metal,  dielectric, 
nonlinear  crystal  etc.,  and  the  atom  is  situated  in  the  vacuum  z  >  0,  with  a 
normal  distance  h  to  the  surface  z  -  0.  Then  the  presence  of  the  substrate 
destroys  the  isotropy  of  the  environment  of  the  atom.  For  instance,  an  emitted 
fluorescent  photon  in  the  -z  direction  can  reflect  at  the  surface  and  travel 
back  into  the  region  z  >  0,  but  an  emitted  photon  in  the  +z  direction  will  never 
hit  the  surface.  Furthermore,  reflection  coefficients  for  media  depend  in 
general  on  the  angle  of  incidence  of  the  radiation,  or  from  the  perspective  of 


■  ■'% 
m 

I 
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the  atom,  these  coefficients  depend  on  the  emission  angle.  The  reflected 
radiation  will  be  experienced  by  the  atom  as  an  external  field,  and  a  stimulated 
transition  can  cause  the  previously  emitted  photon  to  be  absorbed  again.  This 
mechanism  effectively  enhances  the  lifetime  of  the  excited  state,  and  thereby  it 
changes  the  Einstein  coefficient  for  spontaneous  decay.  Due  to  the  loss  of 
spherical  symmetry,  however,  the  lifetime  of  a  particular  state  of  the 

multiplet  is  not  necessarily  independent  of  m^  an3rmore. 

In  this  paper  we  present  a  general  theory  for  atomic  spontaneous  decay  near 
a  surface,  without  regard  to  any  of  the  properties  of  the  medium.  We  shall  only 
use  the  remaining  symmetries  of  the  system.  In  this  fashion  we  can  disentangle 
the  pure  geometrical  features  of  spontaneous  decay  from  the  effects  which  result 
from  particular  optical  properties  of  the  substrate.  We  shall  only  assume  that 
the  medium  is  isotropic,  which  pertains  to  most  practical  situations.  The  first 
symmetry  which  remains  is  the  invariance  of  the  system  for  rotations  around  the 
z-axis.  For  a  rotation  over  an  angle  a,  equation  (1.2)  reduces  to 

-im  a 

Pr(«,0,0)  |jgm^>  -  e  ® 

e.g.,  the  state  |jgrag>  transforms  into  itself  (apart  from  a  phase  factor). 
Therefore,  this  symmetry  does  not  give  any  information  about  the  m^-dependence 
of  a  lifetime.  A  second  symmetry  is  the  invariance  for  reflections  in  a  plane 
through  the  z-axis.  If  we  take  this  plane  (arbitrarily)  as  the  xz-plane,  then  a 
reflection  in  this  plane  is  equivalent  to  the  product  operation  of  a  rotation 
over  K  around  the  y-axis,  followed  by  a  parity  operation  (point  reflection  in 
the  origin) .  Then  we  have 


m  > 
e  e 


PxPj^(0,jr,0)  I  jgnig> 


(1.4) 


with  the  reflection  operator  and  P  the  parity  operator.  With 


e  e 


1  -m 

(-1)  ®  ®  5 
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(1.5) 


we  find  from  equation  (1.2) 


R  j  m  > 
xz  ‘  e  e 


m  > 
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(1.6) 


where  the  sign  depends  on  the  parity  of  Consequently,  the  states 

ljgni^>  and  experience  the  same  electromagnetic  environment,  and  hence 

their  lifetimes  are  identical.  If  we  denote  the  Einstein  coefficient  of  the 

state  |i  m  >  by  A  then  we  must  have 
'  e  e  m 
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A 


-m 
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(1.7) 


as  a  ressult  of  the  reflection  symmetry. 

Apart  from  these  symmetries,  we  have  a  causality  requirement  which  imposes 

restrictions  on  the  values  of  A  If  the  atom  is  far  away  from  the  surface, 

e 

every  Einstein  coefficient  A  must  reduce  to  its  free -space  value  A-. 

m  I 

e 

Therefore  we  have  the  condition 


A  -♦  A^,  for  h  -♦  « 
m  I 
e 


(1.8) 
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To  see  this,  we  recall  that  spontaneous  decay  is  intimately  related  to  the 

emission  of  photons.  If  the  distance  between  the  atom  and  the  surface  is 

larger  than  the  lifetime  l/A^^^  of  the  state  |jg™g>>  times  the  speed  of  light, 

e 

then  there  can  be  no  interference  between  the  emission  of  photons  by  the  atom 

and  reflected  photons  by  the  surface  (which  were  emitted  earlier) ,  because  the 

travel  time  of  a  photon  between  atom  and  surface  exceeds  the  emission  time. 

Therefore,  for  distances  h  larger  than  roughly  c/A^  ,  the  atom  behaves  as  an 

e 

atom  in  empty  space.  With  X  —  the  wavelength  of  the  radiation,  we  find 

that  for 


h  >  , 

^f 


(1.9) 


the  surface  effects  should  disappear.  For  low-lying  atomic  transitions  the 

ratio  <"q/A£  is  of  the  order  of  10^,  and  (1.9)  overestimates  the  value  of  h  by 

many  orders  of  magnitude,  for  most  cases.  For  dielectrics  and  metals  we  know 

[3-5]  that  A  =»  if  h  >  A,  due  to  strong  interferences  of  the  various 
e 

incident  and  reflected  waves.  Only  for  very  special  cases,  like  four-wave 
mixing  crystals,  or  phase  conjugators  [6,7]  the  causality  requirement  (1.9) 


imposes  an  actual  upper  limit  for  h,  at  which  surface  effects  should  disappear. 


Spontaneous  decay  is  a  relaxation  phenomenon,  which  is  brought  about  by  the 
coupling  of  the  atomic  dipole  n  to  the  electromagnetic  field.  The  analysis  of 
spontaneous  decay  starts  with  the  full  equation  of  motion  for  the  density 
operator  p(t)  of  atom  plus  radiation 


(2.1) 


i  ^p(t)  -  (L^  +  +  L^pp(t) 

where  the  Liouvillians  represent  the  atom,  the  radiation  and  the 

interaction,  respectively.  They  are  related  to  the  corresonding  Hamiltonians 
according  to 

L^a  -  ,  i  -  a,  r,  ar  ,  (2.2) 

which  defines  their  action  on  an  arbitrary  Hilbert  space  operator  a.  The 

Liouvillian  L  includes  the  modifications  of  the  radiation  field  due  to  the 
r 

presence  of  the  medium.  In  other  words,  H^  is  the  Hamiltonian  for  the  empty 
half- space  z  >  0  (the  electromagnetic  vacuum)  and  the  material  in  z  <  0.  The 
quantity  of  interest  for  spontaneous  decay  is  the  state  of  the  atom, 
iiiuapective  of  the  state  of  the  radiation  field.  This  reduced  atomic  density 
operator  P^(t)  is  defined  by 


Pa(t)  -  Trj.p(t)  ,  (2.3) 

where  the  trace  runs  over  all  states  of  the  radiation  field. 

It  is  a  standard  procedure  is  reservoir  theory  to  derive  an  equation  of 
motion  for  p^(t).  In  the  compact  Liouville  notation  it  reads  [8,9] 

1  -  '‘■a  •  ^‘■>'’a<'>  ■ 

with  L  the  free  evolution  of  the  atom  (no  coupling  to  the  radiation) ,  and  P  the 
Si 

spontaneous -decay  operator.  Explicitly, 


-  Tr  L  r 
r  ar  Jq 


-l(L  +L  )r  iL  r 

j  a  r  ,  a  ,  -  > 

dr  e  L  e  (a  p  ) 

ar  a^r 


(2.5) 


(3.1) 


H  -  -n’Eih)  . 
ar  c 


with  h  -  he^  the  position  of  the  atom. 

The  eigenstates  of  the  atomic  Hamiltonian  (internal  structure)  are  the 

angular  momentum  states  |j  m  >,  m  -  -j  . j  and  |j  m  >,  m  -  -j  .  ■  •  •  . j  ■ 

eee  e  e  S&S  &  6 

The  eigenvalue  equations  are 


H|lm>-)<w|jm> 
a'^'ce  e'-’ee 


Hllm>-V(w|jra> 
a'-'g  g  g'-^g  g 


(3.2) 


(3.3) 
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for  an  arbitrary  atomic  operator  a  .  Here,  p  is  the  thermal-equilibrium 
density  operator  of  the  radiation  field,  which  will  be  assumed  to  be  the  vacuum 
state  1 0x0 1  ,  defined  as  the  lowest-energy  state.  Notice  that  1 0x0  |  is  not 
necessarily  the  same  as  for  a  radiation  field  in  empty  space  (zero-photon  Fock 
state) ,  because  the  medium  in  z  <  0  will  affect  the  state  of  the  radiation  in 
z  >  0 . 


Ill .  Dioole  Interaction 

Expression  (2.5)  for  the  relaxation  operator  F  holds  quite  generally.  An 
explicit  evaluation  (for  instance  its  matrix  elements)  requires  that  we 
prescribe  the  interaction  Hamiltonian  If  denote  by  E(r)  the  electric 

component  of  the  radiation  field  (although  further  unspecified) ,  then  the 
coupling  Hamiltonian  in  the  dipole  approximation  assumes  the  form 


■vW.V.V .V 


and  the  frequency  separation  between  the  two  doublets  is  >  0. 

terms  of  the  projectors  onto  the  e  and  g  levels 


P  -  )  |j  m  Xj  ml  ,  P  -  )  (j  m  xj  m 
e  /  'Je  e  -'e  e'  g  [_  '-Jg  g  Jg  g 


the  Hamiltonian  can  be  represented  by 


H  -  hw  P  +•  P 
a  '  e  e  g  g 


(3.4) 


(3.5) 


Then  the  evaluation  of  the  exponentials  exp(±iL^r)  in  (2.5)  proceeds  in  two 
steps.  From  L  a  -  K  ^[H  ,a]  it  follows  that 

3  Si 


±iL  r  ±iH  r/Jf  +iH  r/>< 
a  a  "  a 

e  <7  -  e  a  e 


2  2 

and  with  P  -  P  ,  P  -  P  we  find 
e  e’  g  g 


±iH  r/M  ±iw  r  ±iw  r 

e  ^  -  e  ®  P  +  e  ^  p 

e  g 


Combining  (3.6)  and  (3.7)  gives 


(3.6) 


(3.7) 


±tL  r 
a 

e  a 


±i(w  -u)  )r 


(3.8) 


cr~e,g 

/9-e,g 


An  important  property  of  an  atomic  dipole  operator  ^  is  that  it  cannot  have 


matrix  elements  between  states  within  a  single  multiplet.  In  terms  of 
projectors  we  can  then  write 

which  is  essentially  Laporte's  rule  (pg.  260  of  Ref.  1).  On  the  other  hand,  the 
closure  relation  for  the  atomic  wave  functions  is 


P  +  P  -  1  ,  (3.10) 

e  g 


so  that  /i-(P  +  P)/i(P  +  P), 

e  e  g'c'  e  g' ' 


and  with  equation  (3.9)  this  reduces  to 


a.  " 


(3.11) 


Here  we  introduced  the  lowering  (+)  and  raising  (-)  part  of  ^  as 

-  P  pP  .  -  F  fiP  ,  (3.12) 

-  gC  e  ’  e-  & 

and  from  -  ju  we  find 

.  (3.13) 

Now  we  substitute  expression  (3.8)  twice  in  equation  (2.5),  and  we  expand 
the  dipole  moment  /i  and  the  field  at  the  position  of  the  atom,  E(h) ,  in 
Carthesian  components.  After  some  rearrangments  we  obtain  the  representation 


in  terms  of  the  Hilbert  space  operators 


j 


(3.15) 


The  nine  scalar  functions  fj^j(r)  (not  operators)  are  the  field  correlation 
functions 


(O 


Tr^  Ej^(h) 


■iL  T 
r 


(E^(h)p^) 


(3.16) 


which  depend  on  properties  of  the  radiation  field  only  (i.e.  E(h) ,  and  p^) . 
With 


E(h,r) 


iL  r 
r 


E(h) 


(3.17) 


the  field  in  the  interaction  picture,  and  with  -  |0X0|,  we  can  write 
equation  (3.16)  as 


fj^j(0  -  <0|E^(h,r)Ej(h.0)|0> 


(3.18) 


which  clearly  exhibits  that  fj^j(r)  is  the  correlation  function  of  the  electric 
field  at  the  space  point  h. 
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Next,  we  insert  expression  (3.8)  into  the  definition  (3.15)  of  Q^,  and  we 


use  (3.9).  The  r-integral  effectively  amounts  to  a  Fourier -Laplace  transform  of 


fj^j(r)  according  to 


dr  e  f  (r) 


(3.19) 


Working  out  expression  (3.15)  then  gives 


Q.  -  )  ?.  .  (w  )P  /i.P 
/_  ly  o'  g^j  e 


(3.20) 


where  we  have  made  the  approximation 


(3.21) 


That  this  is  as  good  as  exact  for  the  electromagnetic  vacuum  follows  fr  "(m  the 


representation  (3.16)  of  fj^j(r).  The  field  (h)  consists  of  a  creation  and  an 


annihilation  part,  but  since  it  works  on  -  |0X0|,  only  the  creation  part 


contributes.  Therefore,  exp( - iL^r)  (Ej  (h)  |  0x0 1 )  contains  mainly  positive 


frequencies.  It  has  terms  like  exp(-iwr)  with  w  >  0  the  frequency  of  a  photon. 


The  integrand  of  (3.19)  then  has  the  factor  exp(-i(uH^^)r)  for  fj^j(-«^),  and 


exp(-i(w-w^)r)  for  I («^)  .  Oscillations  with  twice  the  optical  frequency,  w  + 


<j^,  will  cause  the  integral  of  exp(-i(w+«^)r)f (r)  over  t  to  vanish  almost 


identically,  as  compared  to  the  same  integral  with  exp(-i(w-w^)r) .  This 


justifies  approximation  (3.21)  for  fields  in  the  vacuum  state  (zero 
temperature) . 


We  then  insert  (3.20)  into  (3.14),  work  out  the  commutator,  use  again  (3.9) 


and  (3.10),  and  drop  nonsecular  terms  [10],  which  finally  yields  for  F 


Fa 


ij 


{f,.(w  )P  u.P  a  +  f .  .  («  )*a  P  /i.P  fi.P 
ij"  o"  e^i  g^j  e  a  ij "  o"  a  e^j  g^i  e 


-  (f,.(a)  )  +  l.,(w  )*)P  fi,?  a  P  n,?  ) 
'  ij"  o'  ji"  o'  '  g'^j  e  a  e'^i  g' 


(3.22) 


We  remark  that  the  only  atomic  property  that  comes  in  the  expression  for  F  is 
the  dipole  operator  n.  On  the  other  hand,  the  functions  fj^j(u)^)  embody  all 
necessary  details  of  the  radiation  field,  the  wave  function  | 0>  of  the  radiation 
field,  the  properties  of  the  substrate  (through  the  modification  of  E(h)),  and 
the  dependence  of  the  spontaneous -emission  operator  on  the  atom-surface 
distance  h  (through  E(h)).  So  far,  we  have  only  used  the  fact  that  the 
radiation  field  is  in  its  lowest  energy  state,  and  therefore  expression  (3.22) 
is  a  veiry  general  representation  for  F  of  an  atom  in  an  empty  part  of  space ,  but 
in  the  vicinity  of  active  boundaries. 


4.  Symmetries 

Although  expression  (3.22)  for  the  spontaneous -emission  operator  of  an  atom 
is  a  great  simplification  as  compared  to  the  general  expression  (2.5)  for  a 
relaxation  operator,  it  still  involves  nine  unknown  field  correlation  functions 
(w^)  .  In  this  section  we  shall  show  that  by  imposing  the  symmetry 
conditions,  as  mentioned  in  the  Introduction,  the  number  of  unknown  parameters 
reduces  from  nine  to  two.  For  an  atom  in  empty  space  it  follows  from  the 
isotropy  of  space  that  all  levels  |jg>“g>  must  have  the  same  Einstein  coefficient 
Aj  for  spontaneous  decay,  as  pointed  out  in  the  Introduction.  The  argument 


relies  on  the  rotational  symmetry  only,  and  is  independent  of  the  mechanism  of 

spontaneous  emission  (e.g.,  the  specific  form  of  H  ).  For  the  remaining 

symmetry  of  an  atom  near  a  surface,  however,  we  only  found  relation  (1.7)  which 

reduces  the  number  of  unknown  Einstein  coefficients  A  ,  2i  +  1,  to  i  +  1  (i 

m  ■’e  -’e  -’e 

e 

integer)  or  +  4  (j^  half- integer) .  Since  we  have  already  an  explicit 
expression  for  F,  equation  (3.22),  we  can  apply  the  symmetry  transformations 
directly  on  that  result.  In  this  fashion  we  can  take  advantage  of  the  knowledge 
of  the  details  of  the  interaction  (-/i*E(h)),  rather  than  working  with  general 
symmetry  arguments  only.  Also  the  fact  that  (3.22)  separates  the  field 
properties  ?j^j(tJ^),  from  the  atomic  contribution,  is  particularly 

convenient.  In  the  Introduction  we  discussed  a  rotation  of  the  atom  in  a  fixed 
environment  (the  electromagnetic  vacutam) ,  which  involves  complicated  rotation 
matrices  (equation  (1.2)).  Since  we  know  the  explicit  occurrence  of  the 
radiation  field  in  the  expression  for  F,  we  can  now  equally  well  rotate  the 
vacuum  and  keep  the  atom  fixed.  Then  symmetry  requires  that  F  for  a  rotated 
vacuum  around  the  z-axis  is  identical  to  the  original  F,  and  the  same  procedure 
applies  for  a  reflected  vacuum  in  the  xz-plane. 

Let  us  first  consider  a  rotation  around  the  z-axis  over  an  angle  a.  Then 
the  unit  vectors  transform  according  to 


e'  -  cosa  e  +  sina  e 
-X  -X  -y  , 


e' 

-y 


-sina  e  +  cosa  e 
-X  -y 


(4.1) 


V 


e' 

-z 


e 

-z 


and  the  field  correlations  with  respect  to  the  rotated  basis  are 


-  h’^<0|(E(h,r).ep(E(h.0).ep|0> 


(4.2) 


Then  synnnetry  invariance  requires 


-  q^Cr)  . 


(4.3) 


for  i  —  x,y,z,  j  -  x,y,z  and  for  every  angle  a,  which  yields  a  set  of  nine 
equations.  For  instance,  with  i  -  x,  j  —  z  we  find 


f  (r)  -  cosa  f  (r)  +  sina  f  (r) 
xz  xz^  yz 


(4.4) 


Since  this  must  hold  for  every  a,  we  can  take  a  -  x,  which  gives  q2(’’)  "  0- 

Then  equation  (4.4)  reduces  to  sina  f  (r)  -  0,  and  if  we  then  take  a  -  x/2  we 

yz 

obtain  f  (r)  -  0.  Working  out  the  nine  equations  consequently  gives  the 
yz 

relations 


f  (r)  -  f  (r)  -  f  (r)  -  f  (r)  -  0  , 

xz  zx  yz  '  zy'  ' 


(4.5) 


f  (r)  -  f  i.r)  , 
xx'  yy 


(4.6) 


f  (r)  -  -f  (r)  , 

xy  y* 


(4.7) 


and  there  is  no  restriction  on 


unit-vectors  transform  as 


e  -  e 

-X  -X 


e'  —  -e 

-y  -y 


(4.8) 


and  the  symmetry  invariance,  equation  (4.3),  gives  immediately 


f  (r)  -  f  (r)  -  0 
xy'’  yx. 


(4.9) 


Therefore,  only  the  three  field  correlation  functions 


can  be  nonzero,  and  the  relation  f  (r)  -  f  (r)  reduces  the  number  of 

XX  yy 


independent  quantities  to  two.  For  an  isotropic  vacuum  we  would  find 


additionally  that  f  (r)  equals  f.,_(r),  but  near  a  surface  there  is  no  universal 

XIX  zz 


relation  between  and  ^  ^  J  double 


summation  in  equation  (3.22)  reduces  to  a  single  summation,  which  is  a  great 


simplification. 


5 .  Evaluation  of  F 


In  equation  (3.22)  the  n^'s  are  operators  in  atomic  Hilbert  space,  and  for 


a  further  evaluation  of  F  we  need  the  matrix  elements  of  with  respect  to  the 


angular  momentum  states.  In  the  i  refers  to  a  Carthesian  component  x,y  or 
z,  but  the  matrix  elements  of  p  are  more  conveniently  expressed  in  spherical 


components  with  respect  to  the  z-axis.  In  terms  of  the  spherical  unit  vectors 


-±1  “  "^^Sx  -  ’  So  “ 


(5.1) 


The  second  symmetry  is  the  invariance  of  the  vacuum  for  a  reflection  in  any  I 
plane  through  the  z-axis.  If  we  take  a  reflection  in  the  xz-plane,  then  the 
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we  can  expand  fx  as 


fi  -  ^  ■ 

r 


with  r 


-1,0,1. 


From 


M  we  find 


(-1)% 


(5.2) 


(5.3) 


Then  the  Wigner-Eckart  theorem  [11]  states  that  the  matrix  elements  of  can  be 
written  as 


<1  m 
•'e  e 


^|j 


m  > 
g  g 


75371 


(5.4) 


and  <j  I I j _ra>  follows  after  complex  conjugation,  in  combination  with  (5.3). 

Now  we  insert  the  expansions  (3.4)  for  the  projectors  and  into 
(3.22),  use  the  properties  of  ?..(w  )  as  found  in  the  previous  section,  and  omit 

Lj  O 

the  small  imaginary  parts  of  ?  (w  )  and  f  (w  )  (the  Lamb  shift),  which  gives 

o  z  z  o 

for  r 


Ta 


A  {d  d'a 
T  T  T  a. 


odd' 
a  T  r 


-  2d‘(7  d  ) 
r  a  r 


(5.5) 


Here  we  introduced  the  ' dipole -allowed  raising  operator' 


”  )  (j„“  I*"  I  j  ra  )  1  j  m  Xi  m  I 

r  g  g  e  e  '•'e  e  -’g  g' 

m  m 

e  g 


(5.6) 
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which  has  the  property  that  it  transforms  a  state  |j  m  >  into  |j  m  >,  but  only 

g  g  ®  e 

if  the  Clebsch-Gordan  coefficient  (i  m  Irlj  m  )  is  nonzero,  e.g.,  if  the 

■^g  g  '-^e  e' 

transition  is  dipole  allowed.  The  coefficients  are  defined  by 


*±i  -  i<Jei‘‘iJg>iVy3vi  . 


(5.7) 


(5.8) 


Equations  (5. 5) -(5. 8)  constitute  the  central  result  of  this  paper.  Expression 
(5.5)  is  the  most  general  form  of  the  spontaneous -decay  operator  for  an  atom 
near  a  surface.  It  is  important  to  realize  that  the  entire  operator  in  curly 
brackets  is  parameter  free,  and  therefore  its  structure  is  completely  determined 
by  the  symmetry  requirements.  Furthermore,  the  two  independent  vacuum- field 
correlation  functions  and  the  matrix  elements  of  the  dipole  operator  only  enter 
the  expression  for  F  through  the  two  parameters  A±1  and  Aq.  These  parameters 
are  proportional  to  ^xx^'^o^  ^zz^*^o^’  independent  of 

the  level  structure  (degeneracies)  of  the  multiplets. 


6 .  Equation  of  motion 

In  this  section  we  consider  the  equation  of  motion  for  an  atom  near  a 
surface,  and  the  significance  of  its  solution  will  be  discussed  in  the  next 
section.  The  solution  of  equation  (2.4)  is 


■i(L  -ir)t 


P^(t)  -  e 


P^(0) 


(6.1) 


>  1 

r.'> 


Then  (6.2)  is  substituted  into  (2.4),  and  we  take  all  possible  matrix  elements 
of  the  equation.  We  obtain 


e  e 


(6.4) 


^<j  m  Ip  (t)lj  m'>  -  )  A  (j  m  Irlj  m  )  ( j  m' Ir  |  j  in' )<j  m  I  p  (t)  |  j  m': 
dt  ■’g  g'  a  ■‘g  g  l_  ^  S  S  ®  e  g  g  e  e  -'e  e'  a  '-'e  e 


m  m'r 
e  e 


(6.5) 


7~<j  m  Ip  (t)|j  m>-  -(^A  +  iw  )<1  m  Ip  (c)|j  m  > 

dt  -^e  e'^^a  ■'g  g  ™e  °  ^  e'  ^  g  g 


(6.6) 


— <j  m  Ip  (t)|i  ra  >  -  -(hA  -  iw  )<i  m  Ip  (t)(j  m  > 
dt  -'g  g'  a  '-^e  e  "'e  °  ^  g  ^  ®  ® 


(6.7) 


which  constitutes  a  set  of  ^(jg  +  jg  1)  equations.  The  solution  is 


-h(A  +A  ,)c 


’  e  a '  a '  '  '  e  e 


(6.8) 


<j  m  Ip  (t)|j  m'>  -  <j  ra  Ip  (0)|j  ra'>  +  ) 

Jg  gif'a  /ijg  g  Jg  gi^3  ijg  g  ^ 


0& 

V"  2A  ni  rn 

'  •  > 

^  ,  m  m 
m  m  r  e  e 
e  e 


X  (j  m  Irlj  ra  ) ( j  ra'lrjj  m')<j  ra  Ip  (0)|j  ra'>  ,  (6.9) 

-* g  g  ' e  e  '’g  g  ’■'e  e  •'e  e'  a  '  e  e 


-(Ha  +iw  )t 
m  o 

<1m|p(t)|jra>-e  ®  <jra|p(0)|jra> 

-'e  e'^^a'  '■^g  g  -’e  e''^a  '  g  g 


(6.10) 


-  ( HA  -  iw  )  t 

<Jmlp(t)|jra>-e  ”e  ^  <jnilp(0)|jra> 


■g  g‘  a'  '  '•'e  e 


'g  g‘"a'  '-^e  e 


(6.11) 


for  any  intlal  state  p  (0) .  In  the  limit  t  -►  ®  all  matrix  elements  approach 

zero,  except  <j  m  |p  (‘«>)|j  m'>,  which  reflects  the  fact  that  all  population  is 
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in  the  ground  state.  Furthermore,  we  notice  that  this  long-time  solution  is  not 


unique,  due  to  the  term  <j  m  \p  (0)|j  m'>  on  the  right-hand  side  of  (6.9).  The 

s  s  ^  s  s 

solution  for  t  ->  ®  depends  on  the  preparation  of  the  system  at  t  -  0. 

From  (3.14)  we  readily  derive 

Tra(r«7a)  -  0  ,  (Fo^)^  -  Fa^  ,  (6.12) 

for  every  a The  first  identity  giiarantees  the  conservation  of  trace  in  the 
time  evolution  of  P^(t:)  ,  and  the  second  relation  implies  that  a  Hermitian 
initial  density  operator  P^(0)  remains  Hermitian  in  its  time  evolution.  Both 
properties  can  easily  be  verified  for  the  solution  (6 . 8) - (6 . 11) . 


7 .  Ej.nstein  coefficients 

In  the  previous  section  we  introduced  the  parameters  A  ,  which  are  the  m  - 

“e 

dependent  Einstein  coefficients  for  spontaneous  decay  of  the  level  |jgni^>  to  any 

ground  state.  This  interpretation  follows  immediately  from  (6.8),  which  becomes 

for  ra  -  m' 
e  e 


The  population  of  m  decays  with  a  lifetime  1/A  ,  and  from  (6.5)  we  can 

e 

determine  to  which  ground  levels  the  population  decays.  We  find  for  m  -  m' 


dt'^^g“gl 


p  (t) I j  m  > 
'Jg  g 


-  )  A  (j  m  Irlj  m  )^<j  m  Ip  (t)li  m  > 
/  r  ■'g  g  '-^e  e  "^e  e'  a  '-'e  e 

m  r 
e 


(7.2) 


showing  that  |j  m  >  loses  its  population  to  |j  m  >  at  a  rate 
®  ®  6  8 


( j gUiglr I j Summed  over  all  possible  (three  at  most)  values  of  this 


gives  A  from  equation  (6.3).  (The  Clebsch-Gordan  coefficient  (j  m  iTjj  m  )  is 

“g  88  ®  ® 


only  nonzero  for  m  +  r  -  m  ,  so  that  the  double  sums  in  (6.3)  and  (7.2) 
8  ® 


effectively  reduce  to  single  sums).  The  splitting  of  A^  in  contributions  from 


different  transitions  is  called  branching,  which  is  illustrated  in  Fig.  1. 


From  (5.2)  and  (5.4)  we  see  that  a  factor  (j  m  lr|j  m  )  comes  from  the 

8  8  e  e 


component  of  the  dipole  operator.  Since  the  lie  in  the  xy-plane  and 

ej^  is  perpendicular  to  that  plane,  we  can  identify 


A 


±1 


.  Ag  -  Ai 


(7.3) 


where  A||(Aj_)  is  brought  about  by  the  parallel  (perpendicular)  component  of 
For  atoms  in  empty  space  we  have  ”  ^zz^'^o^’  therefore  A||  -  Aj^  (see 

(5.7),  (5.8)),  and  they  both  equal  A^  from  (1.1).  With  the  sum  rule  for 
Clebsch-Gordan  coefficients 


)  (jmlrljm)^-! 

L  8  g  ®  ® 

m  r 
8 


we  then  find  from  (6.3)  that  A  -  A^.  for  every  m  .  Conversely,  the  m  - 

e 

dependence  of  A  for  atoms  near  a  surface  comes  from  the  possibility  that 
e 

?^(<4>^)  is  not  necessarily  equal  to  lack  of  spherical 

symmetry. 


With  the  help  of  (7.4)  we  can  also  perform  the  summation  in  (6.3),  which 


gives 


n M lUULPUWVJVIIWW; ^TOrst^'  -;.■  X.' ^.■'.  V^r^.'^^.T■L■*■■.’^^.'■~■.'^.■.■‘' y;  WJ vv  v/v^ 


m 


A|l  +  (j^m^lOl  j^m^)  (Aj.-A||) 


(7.5) 


This  shows  again  that  A^  is  independent  of  whenever  Aj_  equals  Ajj .  From  the 

e 

properties  of  Clebsch-Gordan  coefficients  we  can  easily  prove  that 

2  2 

0|j,  depends  only  on  via  which  implies 


A  -  A 
-m  m 
e  e 


(7.6) 


This  is  the  symmetry  relation  which  follows  from  reflection  invariance,  as 
pointed  out  in  the  Introduction.  Then  we  recall  the  relation 


M  ,1-  >2  ' 

<JgV'l2A>  -  3 


(7.7) 


m  m 
e  g 


Summing  equation  (6.3)  over  ra^  then  gives 


U  \  -  3  *  t 


(7.8) 


m 


Since  there  are  21  +1  values  m  and  three  values  r,  (7.8)  expresses  chat  the 

e  e 

average  of  A  equals  the  average  of  A  Furthermore,  we  like  to  emphasize  that 
e 

the  m  -dependence  of  A  is  mainly  geometrical  and  independent  of  the  details  of 
e  m^ 

the  radiation  field  in  the  vacuum.  This  follows  from  (6.3),  which  has  only  the 
A||  and  Aj_  as  parameters.  The  m^-dependence  comes  from  a  Clebsch-Gordan 
coefficient  and  is  therefore  model  independent. 

Finally,  we  recall  the  causality  condition  which  states  that  for  a  large 


a 


I? 


V 

V 


w", 

iT 

$ 

y. 


If 


atom-surface  separation  the  's  must  reduce  to  their  free-space  value  A^. 

e 

we  introduce  the  parameters 


-  A^/Aj  ,  r  -  -1,0,1  ,  (7.9) 

and  similarly  b|j  and  bj_,  then  (6.3)  can  be  written  as 

A  -  A„  )  b  (j  m  Irh  m  )^  (7.10) 

nig  f  /  r'-Jg  g  '-Je  e' 

m  r 
g 

The  two  dimensionless  parameters  b||  and  bj_  determine  completely  the  spontaneous - 
decay  operator.  They  depend  on  the  details  of  the  radiation  field,  incorporate 
the  h-dependence  of  T,  and  contain  all  necessary  information  about  the 
substrate,  like  its  dielectric  constant  or  reflectivity.  Everything  else  is 
determined  by  symmetry.  By  definition,  the  b^'s  obey 

lim  b  -  1  (7.11) 

h-Ho 


8.  Determination  of  b[^  and  b|| 

So  far  we  have  only  applied  symmetry  considerations,  and  it  appeared  that 
the  spontaneous -decay  process  is  completely  determined  by  the  two  dimensionless 
parameters  bj|  and  b^.  Evaluation  of  these  quantities  requires  additional 
knowledge  about  the  optical  properties  of  the  medium,  or  about  the 
electromagnetic  vacuum  field.  In  this  section  we  illustrate  with  simple 
examples  how  bit  and  bj^  can  be  computed  in  specific  situations. 


8.1.  Direct  method 


From  (7.9)  and  the  definition  of  we  find  the  explicit,  defining 
relations 


6ir«  J<c 

b||  -  - 2 -  f  ) 

"  3  XX '  o ' 


(8.1) 


b,  - 


6ir«  Wc' 
o' 


?  (w  ) 

zz  o 


(8.2) 


in  terms  of  the  field  correlation  functions  from  (3.8).  For  the  calculation  of 


field  E(r)  in  r  -  h,  the  Hamiltonian  of  the  field,  and  the  wave  function  | 0> 
of  the  vacuxim.  It  will  be  obvious  that  in  most  practical  cases  this  information 
is  not  available.  For  the  trivial  case  where  the  medium  is  absent,  we  have 


E(r)  -  2  ^ 

ks 


It  V  ^ks  -k 
o  -  -s 


ik*r 

e  ~  +  H. c 


(8.3) 


with  V  the  quantization  volume,  a^^  the  annihilation  operator  for  photons  in  the 
mode  ks ,  and  t,  a  unit  polarization  vector.  The  Hamiltonian  is 


with  -  ck,  and  | 0>  is  the  zero-photon  Fock  state.  Standard  calculations  [12; 
then  give 


S.  . 

oir  €  Wc 
o' 


do)  u) 

Jo 


3  -  iwr 

e 


(8.5) 


and  we  find  indeed  f  (r)  -  f  (r)  -  f  (r),  and  f..(r)  -  0  for  i  j .  The 
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Fourier -Laplace  transform  at  is 


I.  .(u  )  -  ^ - r  . 

o' 


(8.6) 


where  we  have  dropped  the  imaginary  part.  Then  equations  (8.1)  and  (8.2)  give 


b||  -  bj_  -  1  , 


(8.7) 


as  it  should  for  an  atom  in  empty  space.  For  more  interesting  cases  this  direct 
method  is  probably  intractable,  because  it  requires  the  full  quantized 
electromagnetic  field. 


8 . 2  Indirect  methods 


If  we  would  be  able  to  calculate  a  quantity  which  depends  on  b||  or  bj^  in  an 

independent  way  (thus  without  reference  to  field  correlation  functions),  then  we 

could  possibly  extract  the  value  of  b||  or  bj^  from  this  additional  information. 

Let  us  consider  the  expression  for  the  Einstein  coefficient  A  from  (7.10), 

m 

e 

which  holds  for  every  level  configuration.  For  the  situation  of  a  j  -  0, 

i  -  1  transition,  the  relevant  Clebsch-Gordan  coefficients  are 
e 


(001r|lm)-5 

'  e  rm 


(8.8) 


and  we  see  that  A  -  Acbii  and  A  -  -  A-bi  .  If  we  could  evaluate  the 

m  -±1  f  m  -0  f 

e  e 

Einstein  coefficients  for  the  levels  jl  1>  and  | 1  0>  by  a  different  method,  then 

this  would  give  directly  b|j  and  bj^.  In  turn  this  would  determine  A^  for  every 

other  level  configuration  via  equation  (7.10). 

As  mentioned  in  the  Introduction,  spontaneous  decay  is  accompanied  by  the 

emission  of  fluorescent  photons,  and  in  such  way  that  energy  is  conserved.  From 

equation  (6.4)  with  m^  =•  we  know  that  the  state  loses  population  at  a 

rate  equal  to  A  times  the  population.  A  transition  from  |j  m  >  to  any  of  the 
e 

ground  states  corresponds  to  the  emission  of  a  photon  with  energy  .  and 
therefore  the  energy  gain  of  the  radiation  field  per  unit  of  time  and  at  time  t 
equals 


-  hw  A  <j  m  Ip  (t) Ij  m  >  . 

dt  om  -^ee'a  '-^ee 
e 


(8.9) 


provided  that  only  the  state  |jgni^>  is  populated.  If  we  prepare  the  atom  at 

time  zero  in  the  excited  state  M  ra  >,  then  we  obtain  from  (7.1) 

'-'e  e 


JTT  -A  t 

^  -  Ww  A  e  e 
dt  '  o  m 

e 


(8.10) 


and  the  emitted  energy  after  completion  of  the  decay  is 


w-f 


(8.11) 
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as  it  should  be.  From  (8.10)  we  also  find 

A  (8.12) 

dt '  t-0  '  o  m 

e 

This  implies  that  we  know  A  (and  thereby  bii  and  b^)  ,  as  soon  as  we  obtain  an 

e 

expression  for  the  emission  rate  at  t  -  0,  or,  for  an  atom  with  density  operator 

I i  m  Xl  ml.  In  the  next  section  we  illustrate  this  method  for  a  particular 
'•'e  e  ■'e  e ' 

example . 

9 .  Perfect  conductor 

The  emitted  atomic  fluorescence  is  dipole  radiation,  and  for  an  atom  in 
empty  space  the  electric  and  magnetic  fields  are  well  known  [13].  Now  let  us 
suppose  that  the  plane  z  -  0  is  a  mirror,  e.g,  the  substrate  is  a  perfect 
conductor  like  silver.  Then  the  reflected  field  can  be  found  from  symmetry 
considerations  (method  of  images),  and  the  total  field  in  the  region  z  >  0  is 
the  sum  of  the  free -dipole  field  and  the  reflected  field.  Inside  the  perfect 
conductor  the  field  vanishes  identically.  In  order  to  find  the  radiated  power 
we  look  at  the  field  far  away  from  the  dipole.  If  we  adopt  a  spherical 
coordinate  system  (r,5,^)  with  respect  to  the  z-axis,  then  the  positive 
frequency  parts  of  the  electric  and  magnetic  fields  for  large  r  are  [14] 

2 

E^'^^(r,t)  -  - ^  {m^'^^(t-r/c,cos^)  -  (r*m^'^\ t-r/c , cos^ ) ) r )  ,  (9.1) 

^  4jre  rc 

o 

n  J- 

B^‘*‘^(r,t)  -  2;^^  £  X  m^‘‘'^(t-r/c,cos<J) 


(9.2) 


with  r  -  r/r  as  the  direction  of  propagation.  The  vector  operator  (t,cosd) 
is  given  by 


m^^\t,costf)  -  (t)cos(cij  he  ^cosd)  -  2i/iM^^  (t)sin(w  he  ^costf) 


(9.3) 


where  the  t-dependence  signifies  the  Heisenberg  picture  of  an  operator,  and  ^ | 
and  the  perpendicular  and  parallel  components  of  the  lowering  part  of 

the  dipole  operator  (equation  (3.12))  with  respect  to  the  xy-plane.  The 
operator  m^^^(t,costf)  represents  the  combination  of  an  atomic  dipole  ^  in  r  -  h 
and  its  mirror  image  in  r  -  -h,  including  retardation.  The  argument  w^hc  ^cosp 
equals  half  the  phase  shift  between  the  radiation  emitted  directly  in  the 
direction  r  by  ^  and  the  radiation  that  is  first  reflected  by  the  surface  and 
subsequently  emitted  in  the  r-direction. 

Now  it  is  an  easy  matter  to  compute  the  radiated  power.  The  emitted  energy 
per  unit  of  time  per  unit  solid  angle  0  in  the  direction  r  should  be  defined  as 
[15] 

2  2 

^  <E^'^(r,t)  X  B^‘*'\r,t)  -  B^'^(r,t)  x  (r ,  t)>*r  ,  (9.4) 

^o  _  _  _  _ 

where  a  minus  (-)  field  is  the  Hermitian  conjugate  of  the  corresponding  plus  (+) 
field.  With  (9.1)  and  (9.2)  we  can  write  (9.4)  as 


a  tan 


<E^'\r.t).E^'^^(r,t)> 


(9.5) 


Then  we  recall  that  <...>  in  the  Heisenberg  picture  stands  for  Tr  p  (0)(...), 

^  di 

and  subsequently  we  transform  the  expression  (9.5)  to  the  Schrodinger  picture, 
which  yields 


Sir  €  c 
o 


-  (f*m^  ^ (costf ) ) (f (costf ) ) ) 


(9.6) 


Here,  m^^\costf)  is  the  Schrodinger  representation  of  ( t , cosS )  ,  which 


.(-) 


m 


(  +  )t 


follows  from  (9.3)  with  the  substitution  /i(t)  -►  p,  and  m' 

Result  (9.6)  gives  the  full  angular  distribution  of  the  fluorescence  in  the 
half-space  z  >  0  for  any  state  p  (t-r/c)  of  the  atom.  The  radiated  power  then 

di 


follows  from 


dW 

dt 


dn 


a  tan 


(9.7) 


where  the  integration  extends  over  half  a  unit  sphere  in  z  >  0,  and  around 
r  -  0.  Elementary  integration  gives 


(i} 

dW  _ ^ 


^  Tr^p^(t-r/c){bj.^->.gi^^  +  . 


(9.8) 


3ir«  c 
o 


where  we  introduced  the  parameters 


bj.  - 


cos(2fl)  _  sin(2fl) 

(2^)^  (2^)^ 


1  -  3{ 


(9.9) 


Vi  •.  \r 


T  3  sin( 23)  cos(2fl1 

■  '  ■  2'  (2^,2 


(9,10) 


y3  -  w^h/c 


(9.11) 


We  shall  see  in  due  course  that  bj^  and  bjj  are  indeed  the  two  parameters  which 

determine  the  spontaneous - decay  operator  T,  and  thereby  the  Einstein 

coefficients  A  of  the  atomic  states  M  m  >. 

m  '-^e  e 

The  raising  part  ''  of  the  dipole  operator  /i  can  be  expressed  in  the 
dipole -allowed  raising  operator  d^  from  (5.6)  according  to 


!  72371  I  ■ 


(9.12) 


and  follows  after  a  Hermitian  conjugation.  The  two  terms  t  -  ±1  then  give 

the  parallel  part  of  ^ ,  and  the  r  —  0  term  is  the  perpendicular  component. 
Combining  everything  gives 


Ap  Tr  p  (t-r/c)  )  b  d  d^  , 

dc  'of  a'^a  L 


(9.13) 


where  we  used  (1.1)  for  the  Einstein  coefficient  A^  in  empty  space.  Next,  we 
substitute  the  expressions  for  d  and  d^  and  perform  the  summations  over  the 
magnetic  quantum  nvunbers .  We  then  obtain 


I 

I 

% 

iK' 


Qi* 
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^  w 
dt  "  ^“o 


Af  ^  b^(jgmgLr|j^m^)2<j^m^|p^(t-r/c)|jgm^>  .  (9.14) 


m  m  r 

e  g 


which  involves  only  the  populations  <j  m  \p  (t-r/c)|j  m  >  of  the  states  |j  m  > 

6  6^  06  60 


at  the  retarded  time  t  -  r/c.  If  the  parameters  ,  found  in  this  section, 


would  indeed  be  the  correct  b  's  for  the  Einstein  coefficients,  then  we  could 

r 


use  (7.10),  and  (9.14)  would  reduce  to 


^  -  )<w  )  A  <j  m  Ip  (t-r/c)|j  m 

dt  o  /  m-^ee'a  '-^ee 

A 


(9.15) 


This  is  exactly  equation  (8.9),  summed  over  m^,  and  as  pointed  out  in  section 


8.2,  an  expression  of  the  form  (9.15)  unambiguously  identifies  the  parameters  bj^ 


and  bii.  Therefore,  we  conclude  that  (9.9)  and  (9.10)  are  the  correct  results 


for  bj_  and  bn  for  the  situation  of  a  perfectly- conducting  substrate.  But  then 


we  also  know  the  field-correlation  functions  ?  (w  )  and  ?__(w  ),  according  to 

XX  o  zz  o 


(8.1)  and  (8.2),  which  are  now  found  without  any  knowledge  of  the  quantized 


field  E(r) ,  the  Hamiltonian  H^,  or  the  wave  function  of  the  vacuum,  |0>.  Of 


course,  results  (9.9)  and  (9.10)  can  also  be  obtained  from  the  field-correlation 


functions  in  an  explicitly-quantized  radiation  field  theory  for  an  empty  half¬ 


space  near  a  mirror  [16,17],  or  from  linear  response  theory  [18,19]. 


We  have  studied  the  spontaneous  decay  of  an  atom  which  is  positioned  in  the 


vicinity  of  an  optically- reflective  surface,  but  without  reference  to  any 


specific  characteristics  of  the  medium.  An  expression  for  the  spontaneous -decay 


operator  T  of  a  degenerate  two-level  atom  was  derived  from  symmetry 
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considerations.  It  appears  that  r  only  involves  the  two  dimensionless 

parameters  b||  and  bj_,  which  are  independent  of  any  of  the  atomic  properties,  and 

are  determined  by  the  vacuum-field  correlation  functions  )  and  f  (w  )  at 

the  position  h  of  the  atom.  The  parameters  b||  and  bj^  incorporate  the  h- 

dependence  of  atomic  lifetimes  and  the  details  of  the  characteristics  of  the 

medium  (like  a  dielectric  constant  or  a  nonlinear  susceptibility) .  We  have 

found  that  the  Einstein  coefficients  for  spontaneous  decay  of  the  levels 

depend  on  the  magnetic  quantum  ntimber  m^,  due  to  the  absence  of  spherical 

symmetry.  Nevertheless,  the  m  -dependence  of  A  is  almost  entirely 

e 

geometrical,  as  is  reflected  in  the  existence  of  the  sum  rule  (7.8)  and  the 
symmetry  relation  (7.6). 

Although  we  only  worked  out  the  case  of  a  two-level  atom  near  a  surface,  it 
should  be  obvious  that  the  same  procedure  applies  to  any  configuration  which  has 
a  symmetry  for  rotation  around  an  axis,  and  a  reflection  symmetry  for  a  plane 
through  that  axis.  For  instance,  the  configuration  of  an  atom  in  between  two 
parallel  mirrors,  as  in  a  recent  experiment  [20],  and  a  geometry  with  an  atom 
near  a  spherical  or  ellipsoidal  macroscopic  body  have  these  symmetries.  For 
these  cases  the  expression  for  F  is  exactly  the  same,  as  are  the  properties  of 
the  Einstein  coefficients  and  the  solution  of  the  equation  of  motion. 

Furtheinnore ,  the  spontaneous -decay  operator  for  an  arbitrary  multilevel  atom  in 
a  configuration  with  these  s)nnmetries  can  also  only  depend  on  the  two  parameters 
b||  and  bj^.  This  follows  from  the  fact  that  the  vacuum- field  correlation 
functions  are  independent  of  the  presence  of  the  atom.  It  is  the  S3nnmetry  of 
the  vacuum  which  determines  the  structure  of  the  spontaneous -decay  process. 
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Fig.  1.  An  excited  atomic  state  |jg®g>  decays  to  the  ground  state,  and  the 

lifetime  for  this  process  is  .  The  dipole  selection  rule  m^  -  m^  -1,  ni^  or 

e 

m^  +  1  allows  in  general  three  pathways  for  the  transition  to  the  ground  level, 

as  indicated  in  this  figure.  The  three  branches  contribute  differently  to  . 

—  2  ^ 

Transitions  to  (j  m  ±1>  have  a  rate  constant  AiiCj  m  l+l|j  m^)  ,  whereas  the 
S  ®  II  &  O 

2 

vertical  transition  has  an  inverse  lifetime  equal  to  Aj^Cj^m^l  0|j^m^)  Their 

sum  equals  A  If  one  or  two  of  the  ground  levels  is  not  present,  then  the 
e 

corresponding  Clebsch-Gordan  coefficient  is  zero.  For  instance,  in  a  -  1, 
j  -  0  system,  the  only  ground  level  is  |00>. 
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